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ABSTRACT 

The  thin-walled  torsion  specimen  has  been  analyzed  by  the  finite  element  method  to  deter¬ 
mine  its  usefulness  in  material  testing.  A  particular  geometry,  having  a  short  gauge  section  and 
thick  shoulder  regions,  was  examined  using  material  constants  for  316  stainless  steel.  The  speci¬ 
men  was  analyzed  with  an  eye  toward  determining  its  usefulness  in  approximating  simple  shear 
deformation.  Both  forward  and  single  reverse  behavior  were  examined. 

Classical  isotropic  and  Prager- Ziegler  kinematic  hardening  (using  Jaumann  stress  rates) 
were  used  in  the  analyses.  The  gauge  section  did  not  contain  uniform  stress  components  through¬ 
out,  but  nevertheless,  the  macroscopic  response  of  the  specimen  compared  favorably  with  simple 
shear  results.  Plastic  deformation  was  not  entirely  constrained  to  the  gauge  section  but  it  extended 
some  distance  into  the  shoulder  region.  Consequently,  the  shear  strain  in  the  gauge  section  must 
be  directly  measured  with  an  extensometer  or  a  correction  factor  used  when  converting  the 
applied  machine  rotation  to  gauge  shear  strain.  This  correction  factor  was  determined  from  the 
finite  element  calculations.  In  reverse  twisting,  the  macroscopic  response  of  the  specimen  again 
matched  well  with  simple  shear  results  but  neither  material  law  provided  even  a  qualitative  match 
with  experiment.  The  limitations  of  the  thin-wall  specimen  are  described  and  quantified  in  this 
report. 
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INTRODUCTION 


Determining  the  behavior  of  materials,  particularly  metals,  to  large  deformation  conditions  has  pre¬ 
sented  substantial  experimental  difficulties  from  the  time  of  Trcsca’s  early  metalworking  experiments  to 
the  present.  Primary  difficulties  have  been  in  determining  the  stress  and  strain  accurately  in  an  experi¬ 
mental  test  specimen  that  undergoes  gross  deformation.  A  homogeneous  test  region  is  required  which 
is  free  from  large  stress  or  strain  gradients  and  is  large  enough  to  measure  displacements  from  which 
strain  can  be  inferred.  Tension,  compression,  rolling,  extrusion,  torsion  and  drawing  are  all  types  of  pro¬ 
cedures  that  have  been  applied  to  metals  to  achieve  large  deformation.  Comprehensive  reviews  of  large 

1  *1 

deformation  experiments  arc  available.  ’ 

The  thin-walled  torsion  specimen  has  received  considerable  attention  in  the  literature  in  recent  years 
because  it  offers  the  possibility  of  a  simple  deformation  field  and  homogeneous  stress  and  strain  states. 

A  specimen  which  has  a  wall  which  is  only  a  small  fraction  of  the  radius  of  the  section  will  have  a  near¬ 
ly  uniform  strain  distribution  through  the  wall  thickness.  The  nominal  stress  state  is  given  by  dividing 
the  required  torque  or  axial  (thrust)  load  by  the  cross  sectional  area  and  mean  radius.  The  specimen  is 
known  to  be  unstable  at  large  strains  when  it  is  proportioned  by  conventional  means  having  long  uni¬ 
form  gauge  sections  with  gradual  transition  to  the  gripping  region.  Torsional  buckling  is  a  primary 
mode  of  failure  for  such  a  specimen.  This  problem  can  be  suppressed  by  shortening  the  gauge  length  so 
that  it  is  only  a  fraction  of  the  diameter  of  the  gauge  section.  A  specimen  of  this  type  was  first  pro¬ 
posed  by  Hodieme^  for  use  in  hot  worxing  studies  and  popularized  by  Lindholm  et  al.'*  Figure  1  shows 
a  generic  sketch  of  this  type  of  specimen.  Notice  the  short,  thin-walled  gauge  section  which  quickly 
transitions  to  the  thick-walled  shoulder  region  where  the  specimen  is  gripped.  Table  1  lists  some  of  the 
researchers  who  have  used  this  type  of  specimen  for  large  strain  testing.  Also  shown  are  two  important 
geometric  ratios:  the  gauge  length  divided  by  the  average  diameter  of  the  wall  region,  and  the  average 
diameter  divided  by  the  wall  thickness.  The  larger  these  ratios  are,  the  more  likely  buckling  is  to  occur. 
This  presents  a  trade  off  in  specimen  design  since  a  larger  diameter/wall  thickness  ratio  gives  a  more 
uniform  stress  and  strain  through  the  wall. 

Although  originally  used  to  compare  normalized  flow  stress  behavior  with  tension/compression  re¬ 
sults  (Hecker),  the  torsion  test  offers  the  potential  for  discriminating  material  models  as  illustrated  by 


1 .  Gil  Scvillano,  J.,  Van  Houtte,  P.,  and  Acmoudl,  E.,  Large  Strain  Work  Hardening  and  Textures,  Progress  in  Materi¬ 
als  Science  volume  25,  1981,  pp.  69-412. 

2.  Hecker,  S.S.,  Stout,  M.G.,  and  Hash,  D.T.,  Experiments  on  Plastic  Deformation  at  Finite  Strains,  Plasticity  of  Met¬ 
als  at  Finite  Strain:  Theory,  Experiment  and  Computation,  E.H.  Lee  and  R.L.  Mallett,  eds..  Proceedings  of  Research 
Workshop  held  at  Stanford  University,  June  29  -  July  1, 1981,  pp.  162-205. 

3.  Hodieme,  F.A.,  A  Torsion  Test  for  Use  in  Metalworking  Studies,  Journal  of  the  Institute  of  Metals,  Vol.  91 ,  1962, 
pp.  267-273. 

4.  Lindholm  U.S..  Nagy,  A.,  Johnson,  G.R.,  and  Hocgfeldt.  J.M.,  Large  Strain,  High  Strain  Rate  Testing  of  Copper, 
ASME  Journal  of  Engineering  Materials  and  Technology,  Vol.  102, 1980,  pp.  376-381. 
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Table  1,  Geometry  of  Torsion  Specimens  Used  in  the  Literature 


Researcher  Gauge  Length/Oiameter  Diameter/Wall  Thickness 


Hodieme®  (1962) 

0.40 

5. 

Bailey,  Haas,  Nawab^  (1972) 

0.42 

6. 

Eleiche  and  Campbell^  (1976) 

0.08 

43. 

Senseny,  Duffy,  Hawley^  (1978) 

0.16 

29.-42 

Lindholm  et  al.^  (1980) 

0.23 

16.4 

Lipkin,  Chiesa,  Bammann^  (1987) 

0.24 

15.7 

White,  Bronkhorst,  Anand^  (1990) 

0.31 

26. 

Weerasooriya  and  Swanson^®  (1991) 

0.30 

27. 

Angle 


Figure  1 .  Geometry  of  the  thin-walled  torsion  specimen. 


5.  Bailey,  J.A.,  Haas,  S.L.,  and  Nawab,  K.C.,  Anisotropy  in  Plastic  Torsion,  Journal  of  Basic  Engineering,  March 
1972,  pp.  231-237. 

6.  Eleiche,  A.M.,  and  Campbell,  J.D.,  Strain-rate  Effects  During  Reverse  Torsional  Shear,  Experimental  Mechanics, 
Vol.  16, 1976,  pp.  281-290. 

7.  Senseny,  P.E.,  Duffy,  J.,  and  Hawley,  R.H.,  Experiments  on  Strain  Rate  History  and  Temperature  Effects  During 
the  Plastic  Dfformation  of  Close-Packed  Metals.  ASME  Journal  of  Applied  Mechanics,  Vol.  45,  1978,  pp.  60-66. 

8.  Lipkin,  J.,  Chiesa,  M.L.,  and  Bammann,  D.J.,  Thermal  Softening  of 304L  Stainle.  s  Steel:  Experimental  Results  and 
Numerical  Simulations,  Proceedings  of  IMPACT  ‘87,  Bremen,  FRG,  May  1987. 

9.  While,  C.S.,  Bronkhorst,  C.A.,  and  Anand,  L.,  An  Improved  Isotropic-Kinematic  Hardening  Model  for  Moderate 
Deformation  Metal  Plasticity,  Mechanics  of  Materials,  Vol.  10, 1990,  pp.  127-147. 

10.  Weerasooriya,  T.,  and  Swanson,  R.A.,  Experimental  Evaluation  of  the  Taylor-Type  Polycrystal  Model  for  the 
Finite  Defromation  of  an  FCC  Metal  (OFHC  Copper),  U.S.  Army  Materials  Technology  Laboratory,  MTLTR  9 1-20, 
1991. 
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the  simple  shear  predictions  of  Nagtegaal  and  deJong.'  *  In  the  past  decade  there  has  been  intensive 
work  aimed  at  describing  material  behavior  in  large  simple  shear  and  even  formulating  constitutive 
models  with  the  specific  goal  of  describing  simple  shear  in  a  more  intuitively  acceptable  manner.  In 
light  of  these  results  the  torsion  test  has  a  new  importance  insofar  as  it  provides  an  approximation  to 
simple  shear  deformation. 

The  kinematics  of  the  Hnite  strain  tension-torsion  of  a  thin-walled  tube  have  been  examined  in  de- 
tail  by  McMeeking.  For  a  umformly  deforming  tube,  the  stretching  tensor  (symmetric  part  of  the  ve¬ 
locity  gradient)  can  be  written  in  curvilinear  coordinates  as: 


i 

t 


0  0 


D  = 


0  - 

r 


(1) 


Here  the  1,  2,  and  3  directions  are  in  the  radial,  hoop  and  axial  directions,  respectively.  The  torsion 
test  provides  an  approximation  to  simple  shear  only  as  well  as  it  restrains  the  changes  in  wall  thickness, 
t,  mean  radius,  r,  and  axial  strain  ,  e,  to  be  zero.  The  geometry  shown  in  Figure  1  attempts  to  enforce 
the  gauge  section  to  be  free  from  normal  stretches  by  the  presence  of  the  large  shoulders  near  the  gauge 
region.  The  massive  shoulders  prevent  large  radius  change.  The  grips  restrain  axial  motion  to  within 
the  stiffness  of  the  testing  machine  frame  and  the  thick  walls  of  the  shoulder  region  transmit  this  axial 
stiffness  to  the  gauge  region.  The  effectiveness  of  the  geometry  of  the  torsion  specimen  to  restrain  the 
radial,  hoop  and  axial  straining  in  the  testing  of  metals  is  one  of  the  primary  features  to  be  determined  in 
evaluating  the  torsion  test. 

The  experimental  programs  referred  to  in  Table  1  all  use  some  form  of  the  thin-walled  specimen  but 
not  in  a  careful  way  for  determining  multiaxial  stress  response.  The  largest  use  is  to  determine  just  the 
shear  stress  response  for  finite  twisting.  This  can  be  meaningful  or  not  depending  upon  the  restraint  in 
the  axial  direction.  For  the  torsion  specimen  to  be  able  to  approximate  simple  shear  it  must  be  able  to 
ensure  axial  straining  in  the  gauge  section  is  very  small.  This  requires  both  stiff  testing  machines  and 
stiff  grips.  Measurement  of  the  induced  axial  force  required  to  prevent  axial  straining  is  desired  in  this 
type  of  test. 

1 1 .  Nagtegaal,  J.C.,  and  de  Jong,  J.E.,  Some  Aspects  of  Non-lsotropic  Workhardening  in  Finite  Strain  Plasticity,  Plas¬ 
ticity  of  Metals  at  Finite  Strain:  Theory,  Experiment  and  Computation,  E.H.  Lee  and  R.L.  Mallctt,  cds..  Proceedings 
of  Research  Workshop  held  at  Stanford  University,  June  29  -  July  1,  1981,  pp.  55-101. 

12.  McMeeking,  R.M.,  The  Finite  Strain  Tension  Torsion  Test  of  a  Thin-Walled  Tube  of  Elastic-Plastic  Material, 
Intemalional  Journal  of  Solids  and  Structures,  Vol.  18,  1982,  pp.  199-204. 
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The  first  experimental  results  which  presented  measured  axial  stresses  from  thin-walled  tubular 
tests  were  from  White  and  Anand  in  1986  for  1 100  aluminum  (see  Aifantis*^)  and  Lipkin  et  al.*  for 
304L  stainless  steel.  Since  these  first  tests  only  a  small  number  of  others  have  been  reported.  White  et 
al.^  give  the  shear  and  axial  normal  stress  histories  for  monotonic  twisting  of  316  stainless  steel,  1 1(X) 
alummum,  and  three  carbon  steels  1020,  1045,  and  1095.  Weerasooriya  and  Swanson^®  presented  simi¬ 
lar  results  for  the  finite  twisting  of  polycrystalline  copper.  In  all  of  the  above  reported  tests  the  speci¬ 
mens  were  twisted  to  a  final  shear  strain  between  0.9  and  5.  The  shear  stresses  invariably  increased 
monotonically  with  deformation.  The  axial  normal  stresses  were  compressive  and  increased  monotoni- 
cally  in  magnitude  for  all  of  the  materials  except  for  copper.  It  had  a  compressive  maximum  then  re¬ 
versed  and  even  became  tensile. 

Lipkin  and  Lowe^'^  presented  the  interesting  result  for  the  reverse  twisting  of  a  specimen  of  304L 
stainless  steel.  The  shear  stress  exhibits  the  same  behavior  as  the  tension/compression  reverse  straining 
behavior  would  be  expected  to  show  (with  the  anomolous  inflection  as  the  stress  passes  through  zero, 
probably  due  to  the  method  of  gripping  the  specimen  with  pins).  The  axial  stress  magnitude  quickly  re¬ 
covered  toward  zero  as  reverse  twisting  commenced  but  then  passed  through  a  nonzero  minimum  and 
gradually  increased  back  toward  the  value  it  had  prior  to  unloading.  This  behavior  is  very  interesting. 

It  is  similar  to  the  multiaxial  rate  betting  seen  in  a  variety  of  studies*^  where,  for  instance,  a  tube  is  twist¬ 
ed  back  and  forth  in  a  cyclical  torsional  manner.  A  transient  in  the  length  change  is  seen  where  the  spec¬ 
imen  initially  shortens  but  then  lengthens  during  each  half-cycle.  Lowe  and  Lipkin'^  have  used  a  poly¬ 
crystal  plasticity  model  to  examine  the  simple  approximation  to  their  reverse  experiment.  They  showed 
correct  prediction  of  some  of  the  qualitative  features  of  the  axial  strain  transient.  They  noted  that  these 
results  were  sensitive  to  the  amount  of  axial  restraint  applied  to  the  calculation.  The  axial  stiffness  of 
the  testing  machine  and  shoulder  region  could  have  an  important  bearing  on  the  material  response. 

PREVIOUS  ANALYSES 

The  thin-walled  torsion  specimen  has  been  used  largely  without  investigation,  either  experimentally 
or  numerically  as  to  its  effectiveness  in  simulating  simple  shear.  The  two  exceptions  are  Johnson*^  and 
Lipkin,  Chiesa  and  Bammann.®  Johnson  conducted  a  two-dimensional  analysis  of  the  test  specimen 
using  the  EPIC-2  lagrangian,  finite  element  hydrocode.  The  analysis  included  both  heat  conduction/gen¬ 
eration  and  strain  rate  effects.  These  results  present  no  information  about  the  stress  and  strain  distribu¬ 
tions  in  the  test  specimen  but  rather  were  concerned  with  matching  a  viscoplastic  constitutive  law  to 


13.  Aifanlis,  E.C.,  The  Physics  of  Plastic  Deformation,  International  Journal  of  Plasticity,  Vol.  3,  1987,  pp.  211-247. 

14.  Lipkin,  J.,  and  Lowe,  T.C.,  Proceedings  of  Plasticity  ‘89,  The  Second  International  Symposium  on  Plasticity  and 
its  Current  Applications,  August  1989,  Tsu,  Japan. 

15.  Freudenthal,  A.M.,  and  Ronay,  M.,  Accumulation  of  Second-Order  Strain  in  Workhardening  Media,  lUTAM 
Symposium,  June  1966,  Springer- Verlag,  New  York. 

16.  Lowe,  T.C.,  and  Lipkin,  J.,  Analysis  of  Axial  Deformation  Response  During  Reverse  Shear,  submitted  to  the  Jour¬ 
nal  of  the  Mechanics  and  Physics  of  Solids,  1991. 

17.  Johnson,  G.R.,  Dynamic  Analysis  of  a  Torsion  Test  Specimen  Including  Heat  Conduction  and  Plastic  Flow, 
ASME  Journal  of  Engineering  Materials  and  Technology,  Vol.  103,  1981,  pp.  201-206. 
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macroscopic  experimental  results  (Lindholm  et  al.**).  This  analysis  does  not  address  the  question  of 
how  well  the  thin-walled  specimen  simulates  simple  shear.  It  is  of  note  that  Johnson  did  observe  that 
not  all  of  the  plastic  deformation  occurs  in  the  thin  test  section.  For  very  large  shear  strains  he  noted 
that  only  87  percent  of  the  specimen  rotation  occurs  in  the  test  section. 

Lipkin  et  al.®  presented  more  detailed  numerical  results  of  their  test  specimen  using  the  DYNA3D 
finite  element  hydrocode.  The  analysis  was  conducted  for  the  express  purpose  of  comparing  experimen¬ 
tal  results  of  the  authors  with  their  particular  constitutive  model.  As  such,  no  comparison  was  made 
with  better  understood,  more  classical  material  models.  This  limits  the  usefulness  of  their  analysis  for 
general  evaluation  of  the  specimen  but  a  number  of  excellent  observations  were  made.  The  gauge 
length  shows  a  slight  lengthening  during  torsion  even  for  perfectly  fixed  specimen  ends.  This  lengthen¬ 
ing  was  both  measured  experimentally  and  simulated  numerically.  The  shear  strain  at  an  element  in  the 
gauge  section  was  shown  to  differ  from  the  average  shear  strain  calculated  from  the  twist  of  the  grips. 
The  shear  strain  was  quite  uniform  in  the  gauge  section  below  shear  strains  of  about  200%  but  as  twist¬ 
ing  was  continued  above  this  level  the  shear  strain  distribution  varied  considerably.  This  variation  was 
in  such  a  way  as  to  suggest  axial  buckling  was  occurring. 

Q 

There  are  two  main  limitations  to  Lipkin  et  al.  analysis.  First,  classical  plasticity  laws  were  not 
used  to  allow  easy  understanding  of  whether  the  macroscopic  numerical  results  yielded  results  consis¬ 
tent  with  the  predictions  of  simple  shear.  Second,  the  finite  element  mesh  that  was  used  was  very  coarse 
containing  only  3  elements  through  the  thickness.  The  elements  they  used  were  eight  node  brick  type  el¬ 
ements  with  reduced  integration.  This  left  only  three  material  integration  points  through  the  thickness 
of  the  specimen  wall.  Also,  no  mesh  convergence  comparisons  were  reported. 

The  results  of  Lipkin  et  al.*  have  shed  new  light  on  the  behavior  of  the  thin-walled  torsion  specimen 
but  have  not  provided  a  detailed,  quantitative  assessment  as  to  the  limitations  and  usefulness  of  the  spec¬ 
imen.  The  current  analysis  addresses  these  questions. 

FINITE  ELEMENT  MODEL 

The  specimen  geometry  that  was  simulated  corresponds  to  the  316  stainless  steel  specimen  reported 
in  White  et  al.^  (See  Figure  1).  The  gauge  length  was  5.9  mm.  The  inside  diameter  was  19.05  mm. 

The  outside  diameter  of  the  gauge  region  was  20.52  mm.  The  outside  diameter  of  the  shoulder  region 
was  38. 1  mm  and  the  shoulder  angle  was  30  degrees. 

The  finite  element  analysis  was  conducted  with  the  ABAQUS  finite  element  program.  All  of  the 
elements  were  eight  node,  linear  displacement  bricks  having  a  full  eight  material  integration  points 
(C3D8  elements  in  ABAQUS). 


18.  ABAQUS  Users  Manual,  1989,  Hibbiu  Karls.son  and  Sorenson,  Inc.,  Providence,  Rl. 
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The  mesh,  shown  in  Figure  2,  discretized  the  specimen  into  just  one  circumferential  slice  of  either 
one  or  five  degree  extent.  A  fine  mesh  within  that  slice  was  used.  Each  slice  had  9  elements  through 
the  wall  thickness  and  15  along  the  axial  length  of  the  one-half  gauge  length.  Again,  the  nodes  on  the 
midplane  of  the  gauge  region  were  restrained  against  both  axial  and  circumferential  motion  but  radial 
motion  was  allowed.  The  nodes  along  the  top  face  of  the  shoulder  region  were  restrained  against  axial 
motion.  The  outer  node  on  the  top  face  was  constrained  to  move  in  a  circular  arc  simulating  the  applied 
twist.  The  compatibility  enforced  on  this  strip  to  make  it  simulate  an  entire  circumference  was  by  re¬ 
quiring  the  corresponding  nodes  on  the  two  faces  of  the  strip  to  retain  the  prescribed  circumferential  an¬ 
gle  (either  one  or  five  degrees)  between  them.  This  required  displacement  relation  between  two  equiva¬ 
lent  nodes  on  opposite  faces  of  the  slice  are  given  by  Equation  2. 


[v] 

X 

COS0  -sin6  0 

V 

y 

=r 

sin  9  COS0  0 

V 

_  0  0  1_ 

(2) 


Here  0  is  the  prescribed  angle  between  faces  A  and  B  and  is  the  global  displacement  in  the  x 
direction  of  the  jth  node  of  the  B  face. 


A  complete  circumferential  mesh  was  used  for  comparison  with  the  single  slice  model.  It  had  a 
coarse  discretization  with  only  7  elements  through  the  wall  thickness  and  18  elements  around  the  cir¬ 
cumference. 


Figure  2.  Finite  element  mesh  used  to  simulate  the  torsion  specimen. 
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The  material  models  that  were  used  in  the  simulations  were  classical,  large  strain  plasticity  laws: 
isotropic  hardening,  kinematic  hardening,  and  perfect  plasticity.  These  were  chosen  since  their  behav¬ 
ior  under  simple  shearing  deformation  is  well  known  (Nagtegaal  and  deJong^  The  Jaumann  stress 
rate  was  used.  It  is  widely  recognized  that  choice  of  the  stress  rate  can  greatly  affect  axial  stress  in  sim¬ 
ple  shear  (Dafalias,*^  Reed  and  Atluri,^®  Aifantts*^).  Finite  element  analyses  using  all  of  the  proposed 
rates  are  beyond  the  scope  of  this  paper.  For  simple  shearing  type  deformation,  the  Jaumann  rate  with 
classical  isotropic  and  kinematic  hardening  bound  the  behavior  both  of  experimental  results  and  of  most 
of  the  constitutive  laws  in  the  literature.  For  elastic  properties,  a  Young’s  modulus  of  200  GPa  and 
Poisson’s  ratio  of  0.33  were  used.  The  material  had  an  initial  yield  stress  of  250  MPa  and  a  constant 
plastic  hardening  modulus  of  1500  MPa.  Convergence  studies  were  conducted  of  the  time  step  size  and 
nodal  force  tolerance  in  the  acceptance  criterion  for  the  finite  element  studies.  The  automatic  load  step¬ 
ping  was  used  in  full  large  deformation  analysis.  Approximately  100  displacement  increments  were 
used  to  twist  the  specimen  to  a  nominal  engineering  shear  strain  of  unity. 

FINITE  ELEMENT  RESULTS  FOR  FORWARD  TWISTING 

The  majority  of  the  comparisons  of  the  numerical  results  are  shown  in  terms  of  macroscopic  vari¬ 
ables  that  would  be  determined  in  an  experiment.  The  nominal  shear  stress  in  the  gauge  section  was  de¬ 
termined  by  taking  the  torque  from  summing  the  circumferential  reaction  forces  on  the  top  face  of  the 
specimen  and  dividing  by  the  cross  sectional  area  and  average  radius  of  the  gauge  section.  The  axial 
normal  stress  was  determined  in  similar  manner  by  summing  the  axial  reaction  force  on  the  top  face  of 
the  specimen  and  dividing  by  the  gauge  cross  sectional  area.  The  average  strains  in  the  gauge  section 
were  determined  assuming  an  extensometer  could  measure  the  circumferential  and  axial  displacements 
on  the  outside  of  the  gauge  section  at  the  intersection  with  the  transition  region  to  the  shoulder.  The 
nominal  axial  strain  was  determined  by  dividing  the  axial  displacement  by  the  original  gauge  length  and 
the  nominal  average  shear  strain  (engineering)  by  dividing  the  arc  displacement  by  the  original  gauge 
length.  These  numerical  results  then  can  be  compared  with  the  variables  measured  in  an  experimental 
:est. 


The  different  models  described  above  were  compared  for  calculations  using  both  isotropic  and  kine¬ 
matic  hardening.  The  macroscopic  stress-strain  results  and  the  stress  and  strain  contours  within  the 
specimen  were  virtually  identical.  The  models:  full  circumferential,  one  degree  slice  and  five  degree 
slice  all  yielded  results  that  were  almost  indistinguishable.  This  was  taken  as  verification  that  the  mesh 
was  sufficiently  fine  and  the  boundary  conditions  for  the  single  slice  meshes  were  appropriate.  Most  of 
the  succeeding  results  that  will  be  discussed  were  obtained  with  the  single  slice  model  having  five  de¬ 
grees  circumferentially. 


19.  Dafalias,  Y.F.,  A  Missing  Link  in  the  Formulation  and  Numerical  Implementation  of  Finite-Transformation  Elas- 
toplasticity.  Constitutive  Equations:  Macro  and  Computational  Aspects,  William,  KJ.,  ed.,  ASME,  1994,  pp.  25-40.. 

20.  Reed,  K.W.,  and  Atluri,  S.N.,  Constitutive  Modeling  and  Computational  Implementation  in  Finite  Strain  Plastic¬ 
ity,  International  Journal  of  Plasticity,  Vol.  1,  1985.  pp.  63-87. 
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Figure  3.  Comparison  of  finite  element  simulation  with  simple 
shear  for  isotropic  hardening  (a)  and  kinematic  hardening  (b). 

In  order  to  evaluate  how  well  the  thin-walled  torsion  specimen  approximates  simple  shear  a  compar¬ 
ison  is  made  between  the  shear  ind  normal  stress  response  inferred  from  the  finite  element  simulation 
with  that  for  the  same  constitutive  models  integrated  directly  assuming  only  simple  shear  deformation. 
The  results  of  this  comparison  are  shown  in  Figure  3  for  both  isotropic  and  kinematic  hardening.  For 
isotropic  hardening  we  see  an  almost  exact  correlation.  No  normal  stress  develops  and  the  shear  stress 
linearly  increases  with  strain.  For  kinematic  hardening  (using  Jaumann  derivatives)  a  substantial  axial 
normal  stress  is  predicted  both  by  the  finite  element  model  of  the  specimen  and  the  assumed  simple 
shear  deformation.  The  finite  element  model  predicts  an  axial  normal  stress  approximately  20%  larger 
than  produced  by  simple  shf  ar.  Notice  that  the  shear  stress  is  quite  well  correlated  by  the  specimen  and 
simple  shear. 

The  distribution  of  various  Cauchy  stress  components  is  shown  in  Figure  4  for  kinematic  hardening 
at  a  nominal  engineering  shear  strain  in  the  gauge  section  of  91%.  Figure  4a  shows  the  cross  section  of 
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the  near  gauge  length  region  as  well  as  the  two  radial  lines,  A  and  B,  that  were  used  to  plot  the  distribu¬ 
tions  of  stress  components  shown  in  Figure  4b.  The  solid  curves  in  Figure  4b  correspond  to  cross  sec¬ 
tional  line  A.  Nonce  that  the  shear  stress  is  quite  uniform  across  the  wall,  thickness  but  the  hoop  and  axi¬ 
al  normal  stresses  have  a  large  variation.  This  is  a  result  of  the  center  of  the  specimen  pulling  in  due  to 
the  tensile  hoop  stress.  For  simple  shear,  the  hoop  stress  should  be  equal  in  magnitude  and  opposite  in 
sign  to  the  axial  stress.  This  can  not  be  maintained  for  a  long  thin- walled  tubular  specimen.  The  diame¬ 
ter  of  the  center  of  the  gauge  section  decreases  slowly  as  the  specimen  is  deformed.  The  diametral  de¬ 
crease  for  the  configuration  shown  in  Figure  4a  is  0.8%.  This  is  not  enough  of  a  change  to  dramatically 
alter  the  geometry  and  would  not  be  considered  a  buckle.  It  does  create  a  superimposed  bending  type  of 
stress  gradient  through  the  wall  thickness  at  this  location. 

The  dashed  curves  show  that  at  location  B  the  stresses  are  much  more  uniform  across  the  wall  thick¬ 
ness.  The  variations  in  both  the  hoop  and  axial  stresses  along  line  B  are  about  16%.  At  this  location  the 
stress  response  of  simple  shear  is  maintained  quite  well. 

The  reason  for  the  increased  axial  normal  stress  when  compared  to  simple  shear  (Figure  3)  can  be 
traced  to  the  slight  diameter  decrease  away  from  the  shoulders  and  the  fact  that  the  plastic  deformation 
is  not  completely  contained  in  the  gauge  section.  Figure  5  shows  the  contours  of  equivalent  plastic 
strain  at  a  grip  rotation  of  0.3  rad.  Notice  that  the  contour  of  0. 1  %  plastic  strain  extends  into  the  transi¬ 
tion  region  for  more  than  one-half  gauge  length.  The  plastic  deformation  is  not  constrained  to  the  gauge 
section.  This  has  two  effects  on  the  kinematics.  The  first  is  that  the  axial  restraint  is  not  perfect. 


Figure  4.  Variation  of  stress  components  with  radial  position  across  the  thickness  at  two  locations 
in  the  gauge,  (a)  Deformed  specimen  geometry  and  location  for  stress  distributions,  (b)  Variation 
of  stress  components  along  line  A  (solid  curves)  and  line  0  (dashed  curves). 
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Figure  5.  Contours  of  equivalent  plastic  strain  plotted  on  the  original  cross  section. 

Figure  6  shows  the  development  of  average  axial  strain  across  the  gauge  length  with  shear  strain. 
For  the  three  plasticity  laws  considered,  an  extension  of  the  gauge  length  is  observed  during  testing. 

The  magnitude  of  this  strain  is  small  but  it  is  enough  to  increase  the  axial  stress,  especially  for  kinemat¬ 
ic  hardening.  For  strict  simple  shear  conditions  in  the  gauge  region,  the  axial  strain  would  be  zero. 

This  small  axial  strain  is  coupled  with  the  slight  decrease  in  the  radius  at  the  center  of  the  gauge  section. 


Figure  6.  Development  of  axial  strain  in  the  gauge  section  with  deformation. 
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Figure  7.  Ratio  of  rotation  in  the  gauge  section  to  total  applied  rotation,  (a)  Variation  of  rotation 
with  axial  position  along  inside  and  outside  surfaces  of  the  specimen,  (b)  Evolution  of  rotation 
ratio  with  applied  rotation  for  various  hardening  laws. 

Another  effect  of  the  plastic  strain  extending  into  the  transition  shoulder  region  is  on  the  shear 
strain.  In  the  results  presented  above,  the  shear  strain  was  calculated  by  taking  the  twist  calculated 
across  the  gauge  section  and  converting  to  strain.  In  experimental  practice,  the  twist  applied  by  the  test¬ 
ing  machine  has  been  assumed  to  be  entirely  transmitted  into  the  gauge  section  and  hence  it  has  been 
used  to  calculate  the  shear  strain.  Figure  7a  illustrates  the  way  that  the  twist  applied  at  the  grips  is  dis¬ 
tributed  along  the  length  of  the  specimen.  Not  all  of  the  rotation  is  confined  to  the  gauge  section  al¬ 
though  the  rotation  does  vary  quite  uniformly  in  the  gauge  section  indicating  uniform  shear  strain.  The 
curves  for  both  the  inner  surface  and  the  outer  surface  of  the  specimen  are  shown.  They  are  identical  ex¬ 
cept  right  at  the  transition  to  the  shoulder  region.  This  would  give  slightly  different  values  for  rotation- 
al  extensometers  (Wu  and  Xu  )  placed  on  the  inner  versus  outer  surface  at  this  location.  From  the  fi¬ 
nite  element  results,  we  compare  the  twist  measured  across  the  gauge  section  with  the  twist  applied 
across  the  entire  specimen.  In  Figure  7b  the  ratio  of  these  twists  (R  =  (A(j))  gauge  /  (A<|))  total)  is  plot¬ 
ted  against  rotation  for  both  the  inner  and  outer  surfaces.  For  the  geometry  considered  here,  the  ratio  is 
seen  to  be  independent  of  hardening  model.  It  is  also  independent  of  applied  strain  after  an  initial  tran¬ 
sient.  The  rotation  ratio  does  depend  slightly  whether  it  is  measured  along  the  inner  surface  or  the  outer 
surface  of  the  specimen.  For  this  specimen,  about  78%  of  the  twist  that  is  applied  at  the  grips  actually 
goes  into  the  deformation  in  the  gauge  section.  It  is  fortuitous  that  this  ratio  is  independent  of  harden¬ 
ing  model  and  deformation  level.  A  simple  correction  factor  can  be  applied  to  the  experimentally  mea¬ 
sured  twist  in  converting  to  shear  strain  (just  multiply  by  0.78).  This  correction  factor  is  dependent 

2 1 .  Wu,  H.C.,  and  Xu,  Z.,  An  Axial-Torsional  Extensometer  for  Finite  Deformation,  ASME  Journal  of  Engineering 

Materials  and  Technology,  Vol.  1 12, 1990,  pp.  330-335. 
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upon  the  particular  specimen  geometry  and  flow  stress  curve  but  can  be  easily  evaluated  from  finite  ele¬ 
ment  modeling.  Of  course,  experimentally  one  would  like  to  have  a  rotational  extensometer  to  measure 
the  twist  in  the  gauge  section.  This  correction  factor  is  the  next  suitable  approach. 

In  order  to  try  to  limit  the  spread  the  of  plasticity  into  the  shoulder  a  simulation  was  considered 
where  the  shoulder  made  an  abrupt  transition  to  the  gauge  with  a  shoulder  angle  of  zero  degrees.  In  re¬ 
sults  not  shown  here  this  square  cornered  specimen  behaved  in  essentially  the  same  manner  as  the  pre¬ 
vious  geometry.  The  macroscopic  stress  response  was  virtually  identical  to  the  results  for  the  tapered 
shoulder.  The  plastic  strain  was  still  observed  to  extend  the  same  way,  as  previously  observed,  into  the 
shoulder  region.  The  results  do  not  seem  sensitive  to  the  angle  of  this  transition.  Note  that  this  applies 
to  the  simulation  with  a  moderate  amount  of  strain  hardening.  The  level  of  strain  hardening  will  deter¬ 
mine  how  well  contained  the  plastic  flow  is  to  the  gauge  region.  For  a  metal  with  little  strain  hardening, 
like  high  strength  steel,  the  plastic  deformation  should  extend  very  little  distance  into  the  shoulder.  For 
a  metal  which  exhibits  a  large  amount  of  strain  hardening,  like  annealled  copper,  the  entire  shoulder  re¬ 
gion  could  deform  plastically  if  not  designed  correctly. 

For  the  tapered  shoulder  geometry  illustrated  in  Figure  1,  the  gauge  length  was  varied  to  see  what  its 
affect  would  be  on  the  axial  stress/strain  history.  The  geometry  considered  so  far  has  a  gauge  length  to 
wall  thickness  ratio  of  approximately  8  (g/t=8).  Simulations  were  conducted  where  the  gauge  length  ra¬ 
tio  was  varied  from  8  to  1.  In  Figure  8  the  macroscopic  stress/strain  results  are  shown  for  these  simula¬ 
tions.  The  trend  can  be  seen  that  as  the  gauge  length  is  shortened  both  the  shear  stress  and  axial  normal 
stress  decrease  in  magnitude.  For  a  ratio  of  4  we  can  see  that  the  simple  shear  results  are  matched  very 
closely.  This  would  be  the  bptimal  ratio  for  this  design  for  this  material.  Unfortunately,  the  analytical 


Figure  8.  Ehect  of  gauge  length  on  the  stress  response  to  kinematic  hardening. 
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simple  shear  results  are  not  matched  for  shorter  gauge  lengths  than  this  optimal  value.  The  usefulness 
of  the  type  of  analysis  shown  here  is  clearly  seen.  The  finite  element  modeling  allows  us  to  choose  the 
best  parameters  for  the  specimen. 

SUMMARY  OF  FORWARD  LOADING  RESULTS 

In  the  preceding  sections  the  thin-walled  specimen  was  analyzed  for  unidirectional  twisting.  The 
main  results  can  be  summarized  as  follows.  Results  were  presented  for  the  stress  strain  distribution  in 
the  gauge  region  and  the  macroscopic  response  of  the  specimen.  Even  though  the  stresses  were  not  uni¬ 
form  throughout  the  gauge  section  the  macroscopic  response  of  the  specimen  was  close  to  that  which 
would  be  predicted  assuming  simple  shear  deformation.  This  good  correspondance  was  attributed  to  the 
existence  of  a  cross  section  in  the  gauge  region  where  the  stresses  were  uniform  and  the  simple  shear  so¬ 
lution  did  hold.  The  calculations  also  showed  that  not  all  of  the  twist  that  was  applied  to  the  shoulder  re¬ 
gion  of  the  specimen  was  transmitted  to  the  gauge  region.  Plastic  deformation  extended  from  the  gauge 
into  the  transition  region  to  the  shoulder.  A  ^actor  was  defined  and  calculated  which  should  be  multi¬ 
plied  by  the  applied  twist  in  converting  the  experimental  torque/twist  curve  to  strcss/shear  strain.  This 
factor  was  constant  over  all  of  the  larger  deformation  and  independent  of  material  model  within  the  lim¬ 
its  bounded  by  classical  isotropic/kinematic  hardening.  This  allows  the  shear  strain  to  be  calculated  cor¬ 
rectly  when  a  torsional  extensometer  is  not  available. 

REVERSE  TORSION  TESTING 

The  thin-walled  torsion  specimen  is  also  examined  as  a  test  specimen  for  conducting  reverse  strain¬ 
ing  experiments.  Deformation  is  conducted  in  one  direction,  then  reversed  to  bring  the  specimen  back  to 
its  original  configuration. 

The  use  of  the  thin- walled  specimen  for  reverse  shear  was  shown  by  Lipkin  and  Lowe^'^  in  1989  for 
304L  stainless  steel.  They  twisted  their  specimen  to  an  engineering  shear  strain  of  about  160%  then  re¬ 
versed  the  twist  back  to  a  total  strain  approaching  zero.  There  was  an  anomaly  in  their  shear  stress  re¬ 
sponse  during  reverse  because  the  specimen  was  gripped  with  a  single  cross  pin.  A  very  interesting  axi¬ 
al  stress  response  was  seen  nonetheless.  During  reverse,  the  compressive  axial  stress  initially  went  to¬ 
ward  zero  very  rapidly  but  then  recovered  and  reached  nearly  the  same  magnitude  that  it  had  prior  to  re¬ 
verse.  This  behavior  was  not  exjjected  since  it  is  not  predicted  by  the  classical  plasticity  laws  of 
isotropic  or  kinematic  hardening  (with  Jaumann  stress  rate).  Lipkin  and  Lowe*'^  used  a  polycrystal  mod¬ 
el  to  produce  this  type  of  qualitative  behavior  although  the  match  with  experiment  was  not  good.  They 
also  observed  a  slight  amount  of  contraction  at  the  center  of  the  gauge  section  which  was  measured  post 
mortem.  They  attributed  it  as  due  to  buckling  during  the  reverse  straining. 

These  interesting  results  of  Lipkin  and  Lowe  serve  as  a  major  impetus  for  this  current  work.  The 
question  of  whether  the  axial  stress  measured  after  strain  reversal  is  really  indicative  of  material  re¬ 
sponse  to  simple  shear  or  whether  it  is  due  to  residual  stress  in  the  specimen  after  forward  straining 
needs  to  be  investigated.  In  other  words,  is  the  thin-walled  specimen  suitable  for  revese  straining  tests? 
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Also,  it  was  desired  to  see  if  other  materials  showed  the  same  interesting  behavior  seen  in  304L  stain¬ 
less  steel.  This  axial  stress  behavior  is  not  predicted  by  the  classical  plasticity  laws  and  it  was  desired  to 
see  what  elements  of  a  phenomenological  constitutive  theory  could  give  this  type  of  behavior.  The  pro¬ 
gram  that  addressed  these  concerns  is  outlined  below. 

First,  the  use  of  the  thin-walled  torsion  specimen  in  reverse  twist  was  investigated  using  the  finite  el¬ 
ement  model  described  above.  The  objective  was  to  see  whether  the  macroscopic  stress  strain  response 
from  the  model  was  the  same  as  the  governing  equations  would  give  for  simple  shear  or  whether  the  ex¬ 
perimental  response  derives  from  some  special  specimen  response.  Second,  reverse  torsion  experiments 
were  conducted  using  the  thin-walled  specimens  of  316  stainless  steel  and  1 100-0  aluminum.  These  ma¬ 
terials  were  chosen  since  a  detailed  data  base  has  already  been  generated  for  them  under  a  number  of  de¬ 
formation  modes  including  torsion  (see  White,  Bronkhorst  and  Anand^).  Finally,  these  reverse  torsion 
tests  were  modeled  using  a  number  of  the  constitutive  laws  presented  in  the  literature  for  modeling  sim¬ 
ple  shear.  Mainly,  these  included  differing  forms  for  the  stress  rate  applied  to  Prager-Ziegler  kinematic 
hardening. 


FINITE  ELEMENT  RESULTS  FOR  REVERSE  TWISTING 

A  finite  element  study  was  undenaken  to  investigate  the  suitability  of  using  the  thin-walled  torsion 
specimen  in  reverse  shear.  This  study  was  an  extension  of  a  detailed  analysis  of  the  torsion  specimen 
discussed  above.  The  same  mesh  and  procedures  described  above  were  also  used  in  this  study.  A  single 
strip  of  8  node,  three  dimensional  brick  elements  was  employed  along  with  suitable  kinematic  con¬ 
straints  to  simulate  a  full  360  degree  circumferential  mesh.  The  strip  model  had  8  elements  along  the  ra¬ 
dial  direction  and  44  elements  along  the  axial  extent  from  the  specimen  symmetry  plane  to  the  end  of 
the  shoulders.  There  were  8x15  elements  in  the  half  gauge  length. 

The  end  of  the  shoulder  region  was  constrained  such  that  it  could  not  move  in  the  axial  direction. 
This  simulated  a  fixed  grip  end  condition. Recall  that  it  was  shown  that  axial  strain  can  still  be  produced 
in  the  gauge  section  due  to  the  induced  stresses  pulling  the  center  of  the  specimen  toward  the  axis.  This 
effect  is  small  and  has  been  quantified. 

The  material  model  in  this  finite  element  simulation  was  linear  kinematic  hardening  using  the  Jau- 
mann  derivative.  An  initial  yield  stress  of  250  MPa  and  a  constant  plastic  hardening  modulus  of  15(X) 
MPa  were  used.  These  were  chosen  as  an  approximation  to  the  compression  test  of  316  stainless  steel  re¬ 
ported  in  White,  Bronkhorst  and  Anand  although  linear  kinematic  hardening  would  not  give  a  good 
match  to  experiment  for  a  uniaxial  reverse  test.  A  Young’s  modulus  of  2(X)  GPa  and  a  Poisson’s  ratio  of 
0.33  were  used  for  the  elastic  properties. 

The  model  of  the  half  specimen  was  twisted  through  an  angle  of  0.3  radians  then  the  direction  of 
twist  was  reversed  and  the  specimen  was  twisted  back  to  the  original  orientation.  The  reaction  forces 
that  were  required  were  summed  to  give  the  torque  and  axial  force.  These  reactions  are  analogous  to  the 
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values  that  would  be  measured  by  a  load  cell  in  an  experiment.  The  macroscopic  specimen  stress  re¬ 
sponse  was  determined  using  the  thin-wall  approximation  by  dividing  the  force  by  the  initial  cross  sec¬ 
tional  area  and  the  torque  by  the  area  and  mean  radius.  The  average  gauge  shear  strain  was  computed  by 
assuming  that  a  torsional  extensometer  could  measure  the  twist  across  the  gauge  section  and  dividing  by 
the  gauge  length. 

The  macroscopic  stress  response  for  this  model  is  shown  in  Figure  9.  The  applied  twist  produced  a 
maximum  shear  strain  in  the  gauge  section  of  75%.  The  shear  stress  shows  the  type  of  response  that  is 
expected  for  a  stress  reversal  with  kinemadc  hardening.  Notice  that  the  strain  magnitude  is  not  large 
enough  to  show  the  stress  oscillations  that  the  Jaumann  stress  rate  would  produce.  The  normal  stress  re¬ 
sponse  shows  the  development  of  a  compressive  axial  stress  during  forward  twisting  and  the  interesting 
result  that  during  reversal  of  the  twist  the  stress  subsides  at  almost  exactly  the  same  rate  that  it  had  de¬ 
veloped.  For  linear  hardening,  the  normal  stress  just  increases  and  decreases  along  essentially  the  same 
curve.  The  finite  element  results  essentially  fall  on  top  of  predictions  where  the  same  material  law  was 
integrated  using  perfect  simple  shear  deformation.  This  demonstrates  the  usefulness  of  the  torsion  speci¬ 
men  for  simulating  simple  shear  under  these  circumstances.  Immediately  it  is  obvious  that  the  shape  of 
the  stress  response  curves  are  very  different  than  those  observed  in  the  experiment  of  Lipkin  and 
Lowe.  The  use  of  the  normal  stress  during  reverse  shear  might  be  a  good  test  for  discriminating  mate¬ 
rial  models.  This  leads  to  the  need  to  investigate  the  shape  of  the  normal  stress  response  during  reversal 
of  shear  from  both  experimental  and  constitutive  modeling  viewpoints. 


Figure  9.  Finite  element  results  for  reverse  torsion.  Solid  curve 
iis  finite  element,  dashed  curve  is  simple  shear  integration. 
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Average  Shear  Strain 

Figure  10.  Axial  strain  development  during  reverse  twist. 


One  last  observation  can  be  made  from  the  finite  element  analysis.  The  small  amount  of  axial  strain 
that  accumulates  during  forward  twist  does  not  recover  during  reverse  twist.  Figure  10  shows  the  aver¬ 
age  axial  strain  across  the  gauge  section  as  a  function  of  the  average  shear  strain.  Notice  that  the  axial 
strain  remains  fairly  constant  during  the  reverse  twist  portion  of  the  curve. 

REVERSE  TORSION  EXPERIMENTS 

Reverse  torsion  experiments  were  conducted  on  thin-walled  specimens  having  short  gauge  sections 
and  thick  shoulder  regions.  The  specimen  geometries  and  testing  procedures  are  similar  to  those  de¬ 
scribed  in  White  et.  al.^  The  inner  bore  of  the  shoulder  region  is  fitted  with  a  matched  machined  steel 
plug  and  the  outer  surface  of  the  shoulder  is  gripped  with  an  hydaulic  collet  type  grip.  This  firmly  holds 
the  specimen  and  allows  both  axial  and  torsional  forces  to  be  applied  with  no  backlash.  The  grips  are 
controlled  so  that  they  experience  no  axial  motion  during  the  twisting  to  within  the  stiffness  of  the  test 
frame.  A  constant  twist  rate  is  applied  in  one  direction  up  until  maximum  rotation  is  reached,  then  the 
twist  rate  is  reversed  and  the  specimen  twisted  back  to  its  original  orientation.  The  twist  rate  was  chosen 
so  that  the  nominal  equivalent  strain  rate  in  the  gauge  section  was  0.001  mm/mm/sec.  No  separate  mea¬ 
surement  was  made  of  the  strain  in  the  gauge  section  apart  of  what  could  be  inferred  from  thegrip  rota¬ 
tion  and  corrected  for  with  the  finite  element  correction.  The  stress  and  strain  results  were  reduced  from 
the  raw  data  in  the  common  way  for  thin-walled  specimens  as  done  above. 

The  results  for  the  two  experiments  are  given  in  Figure  1 1.  The  shear  stress  shows  behavior  very 
similar  to  reverse  loading  of  a  uniaxial  specimen:  plastic  flow  in  the  forward  direction,  elastic  unload¬ 
ing,  a  reduced  yield  stress  level  during  reversing,  a  smooth  elastic-plastic  transition,  and  fully  developed 
plastic  flow  in  the  reverse  direction.  The  axial  stress  response  of  316  stainless  steel  in  Figure  1  la  shows 
behavior  similar,  but  not  identical,  to  the  results  of  Lipkin  and  Lowe.*^  During  forward  loading,  the  axi¬ 
al  normal  stress  magnitude  monotonically  increases  to  100  Mpa.  This  agrees  with  the  magnitude 
reached  in  the  torsion  tests  reported  in  White  et.  al.^  During  reverse  loading,  the  stress  magnitude  quick¬ 
ly  reduces  to  zero  and  then  begins  to  increase  compressively  again.  It  reaches  almost  the  same  magni¬ 
tude  that  it  had  prior  to  unloading  by  the  time  that  the  specimen  is  completely  untwisted.  TTie  results  for 
the  1 100-0  aluminum  are  similar  to  the  316  stainless  steel  but  the  shear  stress  shows  a  region  of  lower 
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strain  hardening  during  the  reversal.  This  behavor  has  been  observed  in  uniaxial  reverse  flow  for  FCC 
materials  that  harden  by  subgrain  formation.  This  is  believed  to  be  the  first  result  that  shows  this  behav¬ 
ior  at  strains  this  large.  The  axial  normal  stress  development  for  the  aluminum  is  similar  in  shape  to  that 
for  stainless  steel  except  that  it  shows  a  litde  inidal  peak  during  forward  loading.  This  could  be  due  to 
an  inital  anisotropy  in  the  aluminum.  Both  materials  were  chosen  and  annealed  to  try  to  obtain  as  ran¬ 
dom  an  initial  crystal  orientation  as  possible  but  this  may  not  have  been  perfect.  No  X-ray  measure¬ 
ments  have  yet  been  made  of  the  crystal  structure  either  before  or  after  deformation.  Note  that  the 
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normal  stress  for  the  aluminum  has  been  multiplied  by  a  factor  of  5  so  that  it  can  appear  on  the  same 
graph  as  the  shear  stress.  The  noise  in  the  load  stress  response  can  be  seen  since  the  load  signal  was  so 
small  for  the  aluminum. 

Detailed  measurements  of  the  geometry  of  the  gauge  section  were  not  made  during  or  after  the  test. 
Some  surface  roughening  was  observed  in  the  gauge  section  as  expected  for  plastic  strains  of  these  mag¬ 
nitudes.  This  was  more  pronounced  in  the  aluminum  than  in  the  stainless  steel. 

PREDICTIONS  OF  REVERSE  SHEAR  WITH  VARIOUS  MODELS 

The  finite  element  results  presented  above  show  that  classical  Prager- Ziegler  kinematic  hardening 
using  the  Jaumann  stress  rate  does  not  predict  the  observed  axial  normal  stress  reponse  in  reverse  tor¬ 
sion.  Many  constitutive  formulations  have  been  proposed  during  the  past  decade  which  modifiy  this  con¬ 
stitutive  law  to  remove  the  unwanted  stress  oscillations  and  match  the  axial  stress  seen  in  experiment. 
Among  these  proposals  are;  the  use  of  the  Green-Naghdi  stress  rate,  and  the  formulation  of  plastic 
spin  in  either  phenomenological^’  or  micromechanical^  frameworks.  The  reader  is  referred  to  these 
papers  for  details  of  the  formulations  and  other  references  on  the  topic. 

In  this  section,  calculations  using  Jaumann,  Green-Naghdi,  and  plastic  spin  motivated  stress  rates  ap¬ 
plied  to  kinematic  hardening  in  simple  shear  are  compared  with  the  experimental  results.  Rigid-plastic 
simple  shear  calculations  are  used  instead  of  full  finite  element  simulations  of  the  test  specimen.  These 
were  much  easier  to  calculate  and  modifiy  than  the  FEM  calculations  would  have  been.  It  is  shown 
above  that  simple  shear  is  an  adequate  approximation  to  the  macroscopic  stress  strain  results  derived 
from  these  test  specimens.  The  plastic  hardening  modulus  used  in  these  calculations  is  derived  from  the 
compression  test  of  316  stainless  steel.  Here,  the  modulus  is  written  as  a  function  of  accumulated  plastic 
strain  and  is  fit  to  the  compression  test  data. 

The  first  prediction  of  theory  with  experiment  is  shown  in  Figure  12.  Here  the  Prager- Ziegler  kine¬ 
matic  theory  using  the  Jaumann  stress  rate  is  compared  to  the  reverse  torsion  of  316  stainless  steel.  No¬ 
tice  that  during  forward  twisting  both  the  shear  stress  and  the  axial  normal  stress  are  predicted  quite 
well  by  the  theory.  During  the  reverse  twist  portion  of  the  results  we  see  large  discrepancies  with  the  ex¬ 
periment.  The  shear  stress  does  not  unload  to  zero  before  fully  developed  plastic  flow  is  predicted  in  the 
reverse  direction.  This  is  to  be  expected  from  any  kinematic  hardening  theory  having  the  Prager-Ziegler 
type  evolution  law  for  the  back  stress.  The  axial  normal  stress  is  predicted  to  monotonically  increase  in 
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(a) 


(b) 


(c) 


Shear  Strain 

Figure  13.  Axial  normal  stress  predictions  for  various  stress 
rates.  Experiment  is  solid  curve,  theory  is  dashed. 


All  of  the  results  presented  above  are  for  Prager-Ziegler  kinematic  hardening.  It  is  natural  to  wonder 
what  the  results  would  be  for  an  evanescent-t)rpe  back  stress  evolut’on  equation.  This  type  of  evolution 
law  was  used  in  White  et  al.^  where  the  material  constants  were  chosen  from  reverse  loading  experi¬ 
ments.  It  is  important  to  remember  that  the  additional  evanscent  term  is  used  to  model  the  elastic-plastic 
transition  during  reverse  cycling.  When  realistic  material  constants  are  chosen  in  this  way  it  was  seen 
that  the  axial  normal  stress  does  not  develop  a  large  enough  magnitude  to  model  the  experiments  in  for¬ 
ward  loading.  Some  calculations  have  shown  that  evanescent  hardening  can  produce  a  qualitatively  bet¬ 
ter  shape  for  the  axial  stress  during  reversal  but  the  magnitude  is  an  order  of  magnitude  less  than  the  ex¬ 
periments. 
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CONCLUSIONS 


These  results  illustrate  the  usefulness  of  the  reverse  torsion  test  for  evaluating  material  models. 

The  results  of  the  finite  element  model  shows  that  the  thin-walled  torsion  specimen  gives  the  same  mac¬ 
roscopic  stress  strain  response  in  both  forward  and  reverse  twist  as  simple  shear.  The  twist  which  the 
grips  apply  to  the  specimen  goes  partly  into  deforming  the  gauge  section  and  partly  into  deforming  the 
shoulder  region.  A  correction  factor  has  been  determined  for  use  in  reducing  the  machine  rotation  into 
the  gauge  region  shear  strain.  This  factor  was  found  to  be  independent  of  the  material  models  consid¬ 
ered  here.  The  stress  and  strain  variations  within  the  gauge  section  were  also  quantified.  The  reverse  tor¬ 
sion  experiments  give  results  that  are  qualitatively  similar  to  those  reported  in  Lipkin  and  Lowe.^'^  The 
axial  normal  stress  initially  decreases  to  zero  but  then  increases  back  to  its  value  prior  to  unloading. 

This  behavior  is  not  predicted  by  Prager-Ziegler  kinematic  hardening  using  any  of  the  corotational 
stress  rates  considered.  Further  research  is  warranted  to  investigate  what  phenomenological  theories  can 
reproduce  this  behavior  without  sacrificing  the  ability  to  model  other  phenomena  such  as  uniaxial  cyclic 
loading.  The  thin-walled  torsion  test  provides  a  new  and  useful  tool  for  modeling  the  deformation  of 
metals. 
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